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We find the most general algebraic type N solution with non-vanishing scalar curvature, which 
comprises all type N solutions of new massive gravity in three dimensions. We also give the special 
forms of this solution, which correspond to certain critical values of the topological mass. Finally, 
we show that at the special limit, the null Killing isometry of the spacetime is restored and the 
solution describes AdS pp-waves. 
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I. INTRODUCTION 



Recently, Bergshoeff, Hohm and Townsend formulated 
a new dynamical theory of massive gravity in three di- 
mensions [l[ . The theory is now called new massive grav- 
ity (NMG) and is described by the Einstein-Hilbert (EH) 
action complemented with a particular higher-derivative 
correction term. The latter provides propagating degrees 
of freedom in the theory, thereby curing dynamically bar- 
ren property of general relativity in three dimensions Q ■ 
In contrast to topologically massive gravity (TMG), for- 
mulated a long ago by Deser, Jackiw and Templeton 
[!, HI, the theory of NMG preserves parity. The asso- 
ciated wave equation contains fourth-order derivatives 
of metric perturbations, representing a physical massive 
graviton with two polarization states. In TMG , the ad- 
dition of the higher-derivative Chern-Simons term to the 
EH action violates parity, but it also makes the theory dy- 
namical with a single propagating massive mode. On the 
other hand, there exist some similarities between TMG 
and NMG theories. First of all, the linearized behavior of 
the latter in the Minkowski background is similar to that 
of TMG, resulting in a unitary theory of two propagating 
massive graviton modes as long as one uses the reverse 
sign for the EH term in the action [l], Q . Attempts to 
extend the NMG theory to all higher dimensions have re- 
vealed that only the three-dimensional model is unitary 
in the tree level 

It is unfortunate that unitarity of both TMG and NMG 
theories requires to reverse the usual sign of the EH term 
in the total action. For instance, due to this pecularity 
the mass of the BTZ black holes |8| becomes negative that 
in turn makes unsatisfactory the quantum description in 
context of the AdS/CFT correspondence. In the case 
of TMG, significant progress on this route was achieved 
recently in @. It was shown that at a "crural" point, 
determined by a certain critical value of the topological 
mass, the bulk gravitons disappear and the BTZ black 
holes have nonnegative masses. That is, one obtains a 
unitary chiral quantum theory of gravity with the usual 
sign of the EH term, which is a dual of two-dimensional 
conformal field theory (CFT2) on the boundary. This 
remarkable result has renewed the interest in TMG with 
the hope of finding other stable vacua for a consistent for- 



mulation of quantum gravity in three dimensions [1 CM - 41 3l j . 
In particular, the authors of [l2| performed the Petrov- 
Segre type algebraic classification of exact solutions to 
TMG, showing that almost all existing in the literature 
homogeneous space solutions locally reduce to either type 
D, biaxially squashed AdS.3 solutions [T^4l6j or type N, 
AdS pp-waves solutions [HI Il7l-20|. New Kundt type 
solutions of TMG were found in 13| . 

A similar wave of activity has also appeared in NMG 
[HHH. Tt was found that NMG admits AdS 3 (BTZ), 
warped AdS3 black holes [HHH as wen as AdS-wave solu- 
tions [22j ■ A new class of asymptotically AdS3 black hole 
solutions to NMG with special values of the cosmological 
term has been discussed in [23], [24{ . Some special cases 
of Kundt spacetimes and homogeneous space solutions to 
NMG were also considered in (23,111]) respectively. In the 
quantum context, contrary to the case of TMG, the rec- 
onciliation of incompatibility of bulk/boundary theories, 
in the sense of their unitarity, still remains unsatisfac- 
tory, though a number of critical relations between the 
AdS3 radius and the mass scale do exist as well HI, [H- 
HH . This fact motivates one to look for all possible exact 
solutions of NMG that could provide stable vacua for 
satisfactory quantum aspects of the theory. 

In a recent paper [32*j . we began an exhaustive pro- 
gramme for studying exact solutions to NMG. In particu- 
lar, we found a simple framework that provided mapping 
all known Petrov-Segre types D and N exact solutions 
of TMG into NMG. Meanwhile, it should be emphasized 
that TMG is a strongly constrained theory as it does not 
admit static solutions besides "trivial" Einstein solutions 
(33], |34j. However, NMG is a much richer theory that ad- 
mits solutions which are absent in TMG [H, [32j . 

In this Letter, we continue the programme of [32j and 
find the most general solution that describes all alge- 
braic type N spacetimes of NMG. In Sec. II we briefly 
recall the field equations of cosmological NMG in terms 
of a first-order differential operator (resembling a Dirac 
type operator) acting on the traceless Ricci tensor. In 
Sec. Ill we introduce a triad basis consisting of two null 
and one spacelike vectors and show that for type N ge- 
ometries, the covariant derivatives of these vectors are 
determined only by three scalar functions. In Sec. IV we 
derive the general form of the spacetime metric with a 
single unknown function of two variables, which obeys 



2 



a second-order linear differential equation. In Sec.V we 
obtain the most general algebraic type N solution and 
discuss its some special limits of interest. 



II. THE FIELD EQUATIONS 

The field equations of NMG were obtained in |l| by 
varying the action 



S = 



16ttG 



1 



d 6 xy/=g [R - 2A ? K 



(1) 



with respect to the spacetime metric. Here R — g^ v R^ u 
is the Ricci scalar, A is a cosmological term, m is a mass 
parameter and 



K = R^R!" 



R 2 



(2) 



In a recent paper [32| , it was shown that if one defines 
a first-order differential operator as 



1 



- 2 ^ 



(3) 



where is a symmetric tensor, then its action on the 
traceless Ricci tensor 



yields 



S j_iu R^V 2 9[tuR 5 



(4) 



(5) 



Here C^„ is the symmetric, traceless and covariantly con- 
stant Cotton tensor defined as 



a 



i 



(6) 



With the quantities p and S^, the field equations of 
NMG with a cosmological term can be put in the form of 
the massive Klein- Gordon type equation with curvature- 
squared source term. Thus, one obtains the field equa- 
tions of NMG in the form 

2 - ra 2 ) = , (7) 

where the traceless source term is given by 
R 1 

This equation is also accompanied by the equation 

R 2 

S^S"" + m 2 R — — = 6m 2 A , 

which is an analogue of the trace equation in [l| . For a 
source tensor given by the relation 



(8) 



(9) 



(10) 



where k is a function of the scalar curvature, which is 
fulfilled for algebraic types D and N spacetimes, instead 
of equation ((TJ), we have 



with 



(11) 



(12) 



We recall that in this description the field equations of 
cosmological TMG acquire the Dirac type form 



(13) 



Here n is the mass parameter of TMG. Further details of 
this description of NMG can be found in 32]. We note 
that in © and in what follows, the Levi-Civita tensor 
is given by the relation e MQ( g = \/—gs tia p and we 
use the convention £012 = 1- 

For algebraic types D and N spacetimes such a 
description turns out to be very powerful for finding 
exact solutions to NMG. This has been demonstrated 
in [32j by mapping all known algebraic types D and N 
exact solutions of TMG into NMG as well as presenting 
new examples of such solutions, which are only inherent 
in NMG. In the following, we further demonstrate the 
advantages of this description and find the most general 
type N solution to NMG. 



III. THE TYPE N GEOMETRIES 

We begin with introducing a triad basis of real vectors, 
{Ip, j n n j m /-i} such that 



1^ — 1 , m M m' 1 = 1 



(14) 



with all other contractions vanishing identically. That is, 
we have two null vectors and and one spacelike unit 
vector m M . The latter vector with 

iv a 

provides an orientation of the three-dimensional mani- 
fold. Clearly, we also have 

(16) 



(15) 



7 1/7(7 V \m® 



The spacetime metric written in terms of the basis vec- 
tors has the form 



(17) 



We now recall that for type N spacetimes, in the Petrov- 
Segre classification of three-dimensional spacetimes, the 
canonical form of the traceless Ricci tensor is given by 



S in/ — Ifili/ , 



(18) 



where Z M is a null vector (see for instance, [ill [2a, l3l 
With this expression, equation (J9j) yields 



A = -v 1 



Am 2 



(19) 
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Here, for a convenience in the following, we have intro- 
duced v 2 — —R/6, which defines the case of a negative 
scalar curvature. Equivalently, one can also define it as 
v = v— A, where A is the usual cosmological constant. 
Meanwhile, using equation (TT0|) we have 

I? = ™ 2 + V — , (20) 

instead of ([12|). 

Next, from the contracted Bianchi identity, S 1 '• " = , 
where the semicolon stands for covariant differentiation, 
we find the relation 

ij% + riw = o. (21) 

This enables one to establish the general representation 
for the covariant derivative of the vector in the form 

l K u = aZ / J„+/3Z /t TO I ,+7 {2m^m u - l^n^ + crm^ly , (22) 

where the coefficients of the expansion are functions of 
spacetime and 7 = . Using this expression, with equa- 
tions (fT5j) and (dH) in mind, one can easily calculate the 
action of the operator p on the traceless Ricci tensor 
S^. We find that 

pS^ = (a- 2/3) y„ - 4 7 l ( „m v) . (23) 

Comparing this equation with that given in (|f 3j) . we see 
that for all type N solutions of TMG the function 7 must 
be zero that, along with (l2"T1) . yields 

Z/ = 0, n^=0. (24) 

That is, the null vector forms a congruence of 
expansion-free null geodesies. In other words, all alge- 
braic type N solutions of TMG are Kundt spacetimes 

B 

Next, we shall show that the above statement remains 
true for the case of NMG as well. For this purpose, it is 
convenient to begin by assuming that 7 7^ 0. Then, the 
use of the Lorentz transformation of the basis vectors 

1 „ 

In -> In , ra^ y TTin flfi , + fm^ - - f Z M , 

(25) 

where / is a real function, enables one to set a — 2/3 in 
the new basis. As a consequence, we have 

1^ = al^lv + j3 {l^m v + 2m^l u ) + 7 (2m lL m v - Z M n„) . 

(26) 

The use of this expression, along with the orthogonality 
condition l^m^ — 0, yields the equation 

m M;v = -2/371^,, - 27 ii M m„ + Z M x„ , (27) 

where is a real vector. It is also straightforward to 
show that 

pS liU = -4r/l (ll m v) . (28) 



With equations (|26[) . (|27l) and (|28|) we are able to calcu- 
late the explicit form of p 2 S^ v , which turns out to involve 
the term proportional to (Z M 7 ;M — 7 2 ) (g^y — 3m^m v ). 
Clearly, this term must vanish, as a consequence of the 
field equations in ([lip . Thus, we arrive at the equation 

^7;m-7 2 - 0. (29) 

We now use the fact that for any vector 

x v\{l > ~ f 1 /);^ — RfivX 11 = (s^y + — g^v j x" , (30) 

where we have used equation (jj). Replacing here x v by 
l v and then contracting the result with Z M , we obtain that 

^7;m + 37 2 = 0. (31) 

Combining now equations ([29| and ([31]) , we see that 7 = 
0, that contradicts with our initial assumption 7 7^ 0, i.e. 
we again arrive at equations in (|24[) . Thus, all algebraic 
type N solutions of NMG are Kundt spacetimes. 

Next, it is important to establish for these spacetimes 
the most suitable representation of the covariant deriva- 
tives of the basis vectors. We first note that 

Pd^ = , l^a = , (32) 

where, with (|24[) in mind, the first equation is obtained 
when using expression (I22[) in (|30[) for x„ — l v and con- 
tracting the result with the basis vector mf 1 , whereas the 
second equation follows from the vanishing divergence 
of equation (f23|) . Equations in (|32|) allows us to reduce 
equation (|22|) into the form 

l K v = odplu + P{lti m v — Ivtrin) . (33) 

by means of the Lorentz transformation of the triad given 
by l^ -> fc x/ % , -t fc _1 ' 2 n M , m M m M , provided 
that l^d^k — 0. We see that equation (|33[) involves only 
two functions a and /3. However, the canonical form of 
the traceless Ricci tensor now acquires an extra function 
k, 

Sfiv — klfjiv , (34) 

not violating the conditions in (|24|) . Using the fact 
that expression (|33[) remains invariant under the Lorentz 
transformation ([25]) , one can choose the vector m M to be 
commuting with the vector in the sense of their Lie 
bracket, 

[l,m]=0. (35) 

This in turn allows us to specify the covariant derivative 
of the vector m in the form 

m^y = Tlfj,ly + (3 {I fitly + lyUfj,) + y\fj,m v . (36) 

With equations ([33p and ([3"6"| one can easily write out 
the covariant derivative of as follows 

= -arifily - /3 (n^my + m^ny) - Ttn^ly - x m ^ m v ■ 

(37) 
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We note that the Lie bracket in (|35|) is preserved with re- 
spect to transformation (I25[) , provided that the function 
/ obeys the condition l^dpf = 0. On the other hand, 
using equations (f3"3")l , (f3l*"[) and ([57]) successively in (f5U|). 
after some manipulations, we also find that l^d^x = 0. 
With these two conditions in mind, one can discard the 
function x in P6|) and (|37[) by means of the Lorentz trans- 
formation in (|25[) . Thus, we find that for the type N 
geometries, the covariant derivatives of the basis vectors 
are determined only by three scalar functions. 

IV. THE CONSTRUCTION OF THE METRIC 

It is straightforward to show that the associated de- 
termining equations for the scalar functions a, /3 and r 
are obtained from equation (|30j) by an appropriate using 
expressions (|3"3")l , (jM)) . (|3"o) and ([57]) in it. As a conse- 
quence, we have the set of simple equations 



d p T = -k- APt , (38) 

d v r = d p a + 2aP = n^dpP , (39) 

d p P = -d v a = v 2 - p 2 , (40) 

d v p = , (41) 



which are easily solved. Here we have used the definitions 
I = l^d^ = d v and m = m^dp = d p . From equation (14*01) . 
we immediately see that for v 2 — (3 2 , the function a is 
constant along the null vector I. Thus, the null vector 
determines a null Killing vector that can be seen from 
Q33p , discarding the function a by a Lorentz transforma- 
tion. This case corresponds to a general AdS pp-waves 
solution [22j (see also below, Sec.V). 

For v ^ /3 2 , the most general solutions to the above 
set of equations are given by 

a = {y 2 - (3 2 ) [-v + f3b(u) + c(u)] , (42) 

r=(u 2 -p 2 ) 2 [vb(u)+g(u,p)}, (43) 

where b(u) and c(u) are arbitrary functions of the coor- 
dinate u, the function g(u,p) obeys the equation 

+ {v2 k p2) 2 = Q . ( 44 ) 

and /3 is determined by equations p0|) and (|4T|) . which 
admit the following solutions 

P = vtanh(vp) , (45) 

and 

P = v coth(^p) . (46) 

Here we have used a coordinate transformation to drop 
a redundant function of u. 
Next, using the Lie brackets 

[n,l] = al, [n,m]=2pn + Tl, (47) 



established by means of equations (|33"j) . ((36]) and (|37)) as 
well as equation ([59^1 . we find the following representation 
for the null vector 

n = n^dp = {v 2 - P 2 ) [A8 V + b{u)d p + d u ] , (48) 

where A is given by 

A = \ [v 2 - a(u, p)} - v [Pb(u) + c(u)] , (49) 

and the function a(u, p) is determined by the equation 

d p a = 2{v 2 - P 2 )g . (50) 

Taking once again the derivative of this equation with re- 
spect to p and combining the result with equations (|38[) 
and (144|) , we obtain the second-order linear inhomoge- 
neous differential equation for a(u, p) 

2k 

d 2 a + 2pdpa = -- T -^. (51) 

The associated dual 1-forms 

Ldx^ = — — du , n u dx^ = dv — Adu , 

v z — p z 

mpdx^ = dp — b(u)du , (52) 

define the metric 
2 

ds 2 = —5 TF^du (dv — Adu) + [dp — b(u)du] , (53) 

v A — P z 

that by means of the coordinate transformation 

v -)• v + Pb{u) + c(u) , (54) 

can also be put in the form 

2 1 

ds 2 = dp 2 H — 7: —r dudv H — 5 -r \a(u, p) — v 2 ] du 2 . 

v z — P z v A — p z 

(55) 

In obtaining this expression, we have used the invariance 
of equations (|4~4"|) and (f5CJ|) with respect to the transforma- 
tion a — > a + Pg\{u) + g2{u), that removes two redundant 
functions in the metric. Thus, the most general metric 
for algebraic type N geometries is characterized by a sin- 
gle unknown function a(u,p) governed by equation (|51[) . 
Clearly, this metric does not admit a null Killing vector 
field. 

The remaining step is to find the explicit form of this 
function. In doing so, as it follows from equations (|51l) . 
we first need to know the function k. For this purpose, 
let us calculate the action of the operator p on the tensor 
in ([M]). We find that 

pSp V = -zSp V , (56) 

where 

z = SP + d p liik. (57) 

Comparing this equation with that in (|13p . we see that 
for z = ±p, we arrive at TMG theory and the corre- 
sponding type N solutions recover those found in (liij . 
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These solutions of TMG can be mapped into NMG using 
the prescription given in [32| . 

It is easy to see that the action of the operator p on 
equation ([55]) yields 



p 2 S l _ lL , = (d p z + z 2 ) , 



(58) 



that, with the field equations of NMG given in (fTTj). leads 
to the equation 

d p z = p 2 - z 2 . (59) 
The nontrival solutions of this equation are given by 

z = ptanh[pp + h(u)] , (60) 

and 

z = p coth[pp + h(u)] , (61) 

where we keep an arbitrary function h(u) as one can not 
gauge out it simultaneously with that entering into the 
solutions for fi (see equations (|45l and (|46l ). 

Next, combining equations (|57|) and (|59|) and taking 
into account equation (|40p . we find that 



k=--{v 2 -fi 2 ) 

2 \n~ — z 



2 fl2\l/2 



2 «2^ | ^ 

2 _ .2 



F(u) , (62) 



where -F(u) is an arbitrary function. 



V. THE GENERAL SOLUTION 

Substituting now the expression of k given above into 
equation (fB"Tj) , making use of solutions (f4"5j) , (|46p and (1501) , 
(|61[) . we can solve it for the metric function a(u,p). 

We first consider the generic case p 2 ^ v 2 and begin 
with fi = i^tanh(z^p) . Then, for the function z given in 
equations (|6H)) and (|6T|) . the solutions of equation ([ST|) 
are respectively given by 



coshfpp + h(u)] 
cosh(^p) 



/(u) + taiih( I /p)/ 1 (u) + / 2 ( U ), (63) 



and 



sinhfpp + „, , , / . „ , , „ , s ,„,s 

o= ^ ; ^ /(u +tanh(^)/ 1 (u ) + / 2 (u). (64) 

cosn(i/p) 

Since h(u) is an arbitrary function, these two solutions 
can be "glued" together to give the single solution 



1 



[cosh(pp) Fi(u) + sinh(pp) F 2 (u) 



cosh(^p) 

+ cosh(z/p) + sinh(^p) f 2 (u)] 



(65) 



which involves four arbitrary functions of it. We recall 
that this solution corresponds to the case of negative 
scalar curvature. However, making an analytical con- 
tinuation v — >• iv or taking the limit v — > , we obtain 
the solutions corresponding to the positive or zero values 



of the scalar curvature, respectively. It is also important 
to note that, as it follows from equation (j2"U)) . the quan- 
tity fx 2 can take on both negative and zero values. The 
associated solutions are also recovered by (|6"5")) when per- 
forming an analytical continuation p — > ip or taking the 
limit p — > 0. 

Similarly, for fi = vcoth.{vp) and p 2 ^ v , we have the 
solution 



1 



sinh(^p) 



[cosh(pp) F\ (u) + sinh(pp) F 2 (u) 



■ cosh(^p) fi(u) + sinh(^p) f 2 (u)] 



(66) 



As in the previous case, one can make appropriate an- 
alytical continuations to include the case of positive or 
zero scalar curvature as well as the case of negative or 
zero p 2 . 

It is interesting to note that the most general solution, 
comprising all these solutions, is given by 



ds 2 = dp 2 



2dudv 
v 2 -fi 2 



Z(u,p) 



v 



du 2 , (67) 



where the metric function has the form 
1 



cosh(pp) F\(u) + sinh(pp) F 2 (u) 



+ cosh(zyp) fi (u) + sinh(^p) f 2 (it)] (68) 



Z(u,p) 



and fi is given as in either ([45]) or (|46|) . We recall that v 
is related to the cosmological term A as given in ([!§)) and 
the general metric does not admit a null Killing vector 
field. For F 2 (u) = fi{u) = f 2 (u) = and fi in flU), this 
solution reduces to that obtained in [32[. We see that 
the general solution involves two extra functions /i(it) 
and f 2 (u). In fact, one of these functions can be dis- 
carded by means of coordinate transformations. In order 
to show this, it is convenient to write equation (|68[) in 
the following alternative form 



Z(u,p) 



cosh(pp) Fi(u) + sinh(pp) F 2 {u) 

hi{u) + fih 2 {u) 
v 2 -fi 2 



(69) 



Here h\{u) — > fi{u) and h 2 {u) — > f 2 (u) for fi given in 
(|45j). whereas hi{u) — > f 2 (u) and h 2 (u) — > fi(u) for 
fi given in (|46j) . Passing now to the new coordinates 
V —> vG{u) + H(u), such that d u G — G(u)H(u), and 
du —> du/G(u), it is straightforward to see that one can 
eliminate h\(u) by choosing G(u) and H{u). As a result, 
we have 



Z(u,p) 



cosh(pp) Fi(it) + sinh(pp) F 2 {u) 
+ sinh(2i/p)F 3 (u) . 



(70) 



Thus, the most general solution is characterized by three 
arbitrary functions. In the following, for some future 
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purposes, we keep both functions fi(u) and f2(u). We 
proceed with the special forms of (|rJT)) , which are of in- 
terest as well. 

(i) p 2 = v 2 (or, as it follows from (f20|). m 2 = v 2 /2). In 
this case, taking properly the limit of (|68| or equivalently 
solving equation (|51|) . with the explicit forms of P and z 
given above, we find that 



Z ( u > P) = i 9 oo { cosh (w) [pFi(u) + h(u)} 
\Jv 2 - p 2 

+ smh(pp)[pF 2 (u) + f 2 (u)}}. (71) 



(ii) p 2 = (or, as it follows from ([20]). m 2 = -v 2 /2). 
Then, from equation (|68[) it immediately follows that 



the metric 

ds 2 = dp 2 + 2e 2up dudv + [e up cosh(fip)Fi(u) 

+e vp sinh(pp)F 2 (u) + e 2vp f^u) + f 2 (u)] du 2 , 

(74) 

where, in the contrary with the general case (|67p. both 
functions fi(u) and f 2 {u) can be gauged out by means 
of coordinate transformations. This metric is nothing 
but the AdS pp-waves solution of NMG that was earlier 
found in 

iv) Finally, for Fi(u) — ±F 2 (u) — f(u), the solution 
in (167p recovers that obtained from the TMG case (see 
Refs.[l3l. |33|). The same remains true for solution (fTTj) 
as well. 

Thus, we have the most general solution given in (|67l) . 
which comprises all algebraic type N spacetimes of NMG. 



Z(u,p) = -=== [Fi («) + P F 2 (u) 
\Jv 2 - p 2 

+ cosh(^p) h(u) + sinhM h («)] • (72) 

Again, one can discard the redundant function fi(u) in 
equations (l7Tj) and ([72]) by means of coordinate transfor- 
mations. 

iii) P 2 — >• v 2 . In this case, making the coordinate trans- 
formation v — > (y — p 2 ) v, one can put the metric in (|6T[) 
in the form 

ds 2 = dp 2 + 2dudv — 4/3w dpdu 

+ [Z{u,p)-{v 2 - p 2 )v 2 ]du 2 . (73) 

We see that in the limit p 2 — > v 2 1 the term proportional 
to v 2 disappears. That is, the Killing isometry of the 
spacetime is restored and the vector d v becomes a null 
Killing vector. This in turn means that the resulting 
metric must describe AdS pp-waves. Indeed, redefining 
once again the coordinate v as v — > e 2up v, we arrive at 



VI. CONCLUSION 

A novel description of NMG in three dimensions, given 
in our previous work [32j , turns out to be a very powerful 
tool for finding exact solutions to the theory. In this 
paper, we have further demonstrated the advantages of 
this formalism. Describing three-dimensional algebraic 
type N spacetimes in terms of a triad basis (with two 
null vectors and one spacelike vector), we have shown 
that all such solutions of NMG are Kundt spacetimes. 
Using this property along with freedoms, provided by 
the Lorentz symmetries underlying the system, we have 
found that the covariant derivatives of the basis vectors 
are in general determined only by three scalar functions, 
which obey a chain of simple equations. 

For algebraic type N geometries, we have obtained the 
general form of the spacetime metric with a single un- 
known function of two variables. Remarkably, this func- 
tion is governed by a second-order linear inhomogeneous 
differential equation. Finally, solving the linear differ- 
ential equation, we have found the most general alge- 
braic type N solution with non-vanishing scalar curva- 
ture, which comprises all type N solutions to the theory 
of NMG. We have also considered some special cases of 
interest. 
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